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ABSTRACT 

The 2004 Dec. 27 giant y-ray flare detected from the magnetar SGR 1806-20 
created an expanding radio nebnla which we have monitored with the Anstralia 
Telescope Gompact Array and the Very Large Array. These data indicate that 
there was an increase in the observed flnx ~25 days after the initial flare that 
lasted for ~8 days, which we believe is the result of ambient material swept-up 
and shocked by this radio nebula. For a distance to SGR 1806-20 of 15 kpc, 
using the properties of this rebrightening we infer that the initial blast wave was 
dominated by baryonic material of mass M > 10^"^'^ g. For an initial expansion 
velocity v ~ 0.7c (as derived in an accompanying paper), we infer this material 
had an initial kinetic energy E > 10^^'^ ergs. If this material originated from the 
magnetar itself, it may have emitted a burst of ultra-high energy {E > 1 TeV) 
neutrinos far brighter than that expected from other astrophysical sources. 
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1. Introduction 

The soft gamma repeater (SGR) 1806-20 is believed to be a magnetar — a slowly 
spinning isolated neutron star with an extremely high magnetic held {B ~ 10^^ gauss; 
Duncan & Thompson 1992, Kouveliotou et al. 1998). On 2004 December 27, a giant hare of 
7 -rays was detected from this object (Borkowski et al. 2004), only the third such event. For 
a distance of 15 kpc (Corbel & Eikenberry 2004; McClure-Grifhths & Gaensler 2005 but see 
Cameron et al. 2005), the Dec. 27 hare was roughly a hundred times more luminous than 
the previous two such events (Palmer et al. 2005; Hurley et al. 2005 and references therein). 
Analysis of a Very Large Array (VLA) observation of SGR 1806-20 seven days after the 
hare discovered a bright, transient source, VLA J180839-202439 (Cameron & Kulkarni 2005; 
Gaensler et al. 2005a), which is believed to have been created by the magnetar during the 
hare. This detection triggered a worldwide radio monitoring ehort, whose initial results have 
been presented by Gaensler et al. (2005b) and by Cameron et al. (2005). In particular, it 
has been determined that the radio source was initially expanding with constant velocity 
V ~ 0.7c (assuming a distance of 15 kpc and one-sided expansion) and that, after day 9, its 
hux decayed as a steep power law (Gaensler et al. 2005b; Taylor et al. 2005). 

Here, we present observational evidence for a short-term rebrightening of this radio 
source which we model as the result of material shocked by ejecta from SGR 1806-20.^*^ We 
then ht the observed huxes to this model, deriving estimates for the mass and energy of the 
ejecta, and discuss this model’s implications for the nature of the Dec. 27 burst. 


2. Observations and Results 

As part of a long-term monitoring campaign of VLA J180839-202439, we have observed 
this source every few days with both the Australia Telescope Compact Array (ATCA) and 
with the VLA. Here we focus on observations at 4.8 GHz from day 6 to day 63 after the 
outburst, as listed in Table 1. The ATCA observations used a bandwidth of 128 MHz, 
and SGR 1806-20 was observed for ~20 minutes at this frequency in each observation. 
For each ATCA observation, we calibrated the flux density scale using an observation of 
PKS B1934-638 at the beginning of the run, and calibrated the phase with a short observation 
of PMN J1811-2055 taken approximately every three minutes. To minimize background 
contamination, we only used data from baselines that included the hxed antenna located 
~ 3 km away from the other five antennae in the array. The VLA observations were reduced 


^°The dynamical properties of this model are described by Granot et al. (2005). 
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using the method described by Taylor et ah (2005) in which hnal phase calibration was 
achieved by self-calibrating the SGR 1806-20 data. For both the VLA and the ATCA 
observations, the radio flux density of SGR 1806-20 was measured using by htting the 
visibility data to a source whose position was a free parameter, htting the visibility data to a 
source whose position is hxed at the location of the SGR, and measuring the peak brightness 
in an image made from these visibilities. In general, these three methods yielded consistent 
results, and any differences are rehected in the errors provided in Table 1. 

The resultant light curve is shown in Figure 1. As reported in Gaensler et al. (2005b), at 
day 9 there was a break in the light curve after which the radio hux faded rapidly. Starting 
on day 15, the observed hux from SGR 1806-20 began to deviate signihcantly from a power 
law decay, and on day 25 the hux began to increase for approximately eight days. On day 
33, the observed hux began to decay again but at a slower rate than between days 9 and 
15. In §3, we model this behavior assuming it is a result of the source’s transition from the 
coasting phase to the Sedov-Taylor phase of its evolution. 

3. A Semi-Analytic Model 

In this Section, we present a semi-analytic model for the evolution of the radio source 
created during the Dec. 27 giant hare. We assume a quasi-spherical shell of hlling factor 
and initial mass M expanding supersonically with an initial velocity Vq into a medium 
of mass density p, driving a forward shock into the ambient material. Initially, the newly 
swept up material is accumulated in a thin layer between the shell and the forward shock, 
and the equation of motion of this shell is: 

^ V = ( 1 ) 

where R = R{t) is the radius of the shell, v = v{t) is the expansion velocity of the shell, and 
P = P{t)^ fhe pressure inside the shell, found from energy conservation to be: 

1 \ f 47r \ 

E = -RRvl = - ( M + —fhR^p] + 2nfbR^p. (2) 

This approximation also works well during the Sedov-Taylor phase (Zeldovich & Raizer 
1966), because even at this stage most of the swept-up material is accumulated in a thin 

^^The results presented by Gaensler et al. (2005b) and Taylor et al. (2005) suggest the radio source is 
elongated and moving along the elongation axis - implying a one-sided outflow and requiring a filling factor. 
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layer just downstream of the shock whereas the rest of the volume is hlled by a rarehed, hot 
gas at nearly constant pressure. By eliminating p and introducing dimensionless variables: 
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At r <C 1, the solution to Equation (4) reduces to u = Uo(l~ 0 . 8 r^). At r 1, the solution to 
this equation asymptotically approaches r = (2.5r^)^/^, close to the Sedov-Taylor solution. 


We assume that, at the forward shock, electrons are heated to an energy yorrieC^ = enipV^, 
where e is proportional to the fraction of the energy density behind the shock in relativistic 
electrons. Electrons with Lorentz factor 7 > 7o are assumed to have a power-law energy 
spectrum N{'y) = ^'( 7 / 70 )“^ (we assume 70 > 1, which is fulhlled for e > 5 x 10“^[c/u]^), 
N{'y)d'y is the number of electrons with energy between yiUeC^ and (y-l-dyjmeC^, K = A^(7o), 
and p is the particle distribution index - which observationally is p ~ 2.5 (Gaensler et al. 
2005b), a typical value for shock-accelerated electrons. Additionally, we assume that the 
magnetic energy density just downstream of the shock front is B‘^/8n = (9/8)es/w^, where 
B is the magnetic held strength, and is the ratio of magnetic to internal energy density 
behind the shock. If the number of emitting electrons is ~ {47r/3)fi,R^p/mp, one can estimate 
the emission from the swept-up material as: 




where d is the distance to the source, S,y is the hux density at a frequency z/, and a = 
4.7 X 10“^® in cgs units. Substituting the above quantities into Equation (5), one obtains: 

Su{t) = llei;^(eB,_in_2)°-®^M24nfo^^dr5^zzGHf/(r) mJy, ( 6 ) 


where dis = d/{15 kpc), the ambient number density n is dehned as n = p/rup and n _2 = 
n/(0.01 cm“^), uio = u/(10^° cm s“^), e_i = e/0.1, es-i = e^/O.!, and the dimensionless 
function /(r) may be found from the solution r(r) to Equation (4): 


/(r) = 



(5p-3)/2 


(7) 


the literature, the electron spectrum is conventionally parametrized by the fraction of the accelerated 
electrons, the fraction of the total energy transferred to these electrons, Ce, and the particle distribution 
index p. In order to avoid cumbersome expressions, we introduce e = ee{p — 2)/2^e(p ~ !)• 
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Both r(r) and /(r) can be found from numerical integration of Equation (4), and are shown 
in Fig. 2. During the coasting phase (r <C 1), the luminosity grows as and reaches a 
maximum at r = 0.78, at which point the expansion velocity has only decreased by 22%. 
At r ~ few, the luminosity decreases as this is faster than the decrease during the 
Sedov-Taylor phase because the pressure within the cavity remains small for a long enough 
time and the expansion velocity decreases faster than in the Sedov-Taylor solution where the 
expanding envelope is hlled by the hot gas. During the Sedov phase (r > 10), the luminosity 
decreases as However, the rate of decline after the maximum depends strongly on the 

microphysics of the shock acceleration (Granot et al. 2005). We do not expect signihcant 
emission from a reverse shock in the ejecta since it was previously shocked by a collision 
with a pre-existing shell (Gaensler et al. 2005b; Granot et al. 2005). 

One can then estimate M and E as: 

M = x 10^"^ S; and (8) 

E = ]^Mvl = 2 . 2 fbn- 24 oViQ x 10^"^ ergs (9) 

assuming the emission peaked SOtso days after the explosion. This estimate for the energy 
is strongly dependent on v, whose uncertainty is dominated by errors in the distance, not 
on projection effects. If SGR 1806-20 was at a lower distance (Gameron et al. 2005), these 
estimates of E and M would decrease signihcantly, though recent results by McGlure-Grifhths 
& Gaensler (2005) support d ~ 15 kpc. Additionally, Uio is related to fb- The expansion 
velocity of n ~ 0.7c quoted in the abstract assumes a one-sided expansion, requiring fb < 0.5. 
Using the elongation observed by Taylor et al. (2005), we derive fb ~ 0.1. 

4. Model Fitting 

To test the model in §3 and to use it to independently estimate the initial mass and 
energy of the source, we £t the observed 4.8 GHz flux densities after day 8.8 to:^^ 

Suit) = So (^ Q + 11 A ^GH^fit/tdec) mJy; (10) 

where Sq mJy is the flux density on day 9, S is the index of the power-law decay, and: 

A = el:^(eB_in_2)°-®^M24nio’^®d75^ (11) 

only used data after day 8.8 in this fit because, as reported in Gaensler et al. (2005b), there is a 
break in the light curve at this epoch which cannot be explained by the model presented in §3. 
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as derived from Equation (6). The fit, shown in Figure 1, was performed using a minumum 
algorithm, and the best-£t parameters (reduced = 1.23) are Sq = 52.4 ± 1.3 mJy, 
5 = —3.12 ± 0.11, A = 11.9 ± 0.2, and t^ec = 46.5 ± 1.7 days. This model predicts that at 
t ~ tdec, the source’s expansion velocity should decrease, as indeed reported at this epoch 
by Taylor et al. (2005). The difference between the observed and predicted shape of the 
rebrightening could be due to several factors — e.g. anisotropy in the outflow (Gaensler 
et al. 2005b; Taylor et al. 2005). However, the fit is good enough that we can use A and 
tdec to express the ejected mass in terms of e, es-i, ■n_ 2 , Uio and di^. Rather than eliminate 
one of these variables, we adopt an expression for M which jointly minimize the power-law 
dependences of all five parameters, finding: 


M = 
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( 12 ) 

(13) 


Here M and E are only weakly dependent on the ambient density, n (which is difficult to 
constrain from observations), but are more sensitive to the shock physics of the flow, e and e^. 
The total energetics of Equation (9) suggest that n _2 < 10^. For di^ ~ 1, n _2 ~ 10, /^ ~ 0.1 
and Uio ~ 2.1 (Taylor et al. 2005), the estimated initial mass is M = 2.1el)’'®^eg°iY x 10^^ g. 
While e and es are unknown, we can estimate them from studies of gamma-ray bursts 
(GRBs) and supernova remnants. If the expanding nebula behaves like the relativistic jets 
produced in a GRB, then e ~ 10“^'® — 10“^'® and ~ 10“® — 10“^ (Panaitescu & Kumar 
2002), implying that M ~ 10^^ — 10^^ g. However, if the behavior of the expanding nebula is 
closer to that of a supernova blast-wave, the magnetic held and relativistic electrons will be 
in energy equipartition, ~ e (Bamba et al. 2003), and e ~ 10“^ — 10“^ (Ellison et al. 
2000), implying M ~ 10^® — 10^^ g. Since it is extremely unlikely that e or is larger than 
0.1, we are rather confident that M > 2.1 x 10^^ g. 


It is also possible that the ambient density is considerably different from n ~ 0.1 cm“^. 
Although the nebula initially expanded into a cavity ~ 10^® cm in size (Gaensler et al. 2005b; 
Granot et al. 2005), by day 25 it had already expanded into the surrounding medium. 
If SGR 1806-20 is inside a stellar wind bubble formed by its progenitor (e.g. Gaensler 
et al. 2005c) or nearby massive stars (Gorbel & Eikenberry 2004), n is possibly ~ 10“^ cm“^, 
implying M ~ 10^^ g. However, SGR 1806-20 is embedded in a dust cloud, n could be 
~ 10 cm“^ implying M ~ 10^® g. In either case, the uncertainty in n does not change the 
order of magnitude of the M and E, which are similar to those derived in Equations (8) and 
(9) that depend on the time of the peak in the light curve (fao) but are independent of the 
shock physics. As a result, we conclude that the Dec. 27 flare created a nebula with initial 
mass > 10^^'® g and initial kinetic energy > 10^^'® ergs. 
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5. Discussion 

An inherent assumption in §3 is that most of the energy of the radio source is in the 
form of modestly relativistic or sub-relativistic baryons, as argued in more detail by Granot 
et ah (2005). We postulate that the source of these baryons is the neutron star itself. The 
giant flare is caused by, and accompanied with, the violent restructuring of the magnetic held 
in which some magnetic held lines may, like a slingshot, throw away the matter from the 
surface layers of the star. Although the canonical picture (Thompson & Duncan 1995, 2001) 
assumes that the magnetic stresses excite predominantly horizontal motions of the crust, one 
can imagine that stretching of magnetic held lines initially buried in the crust may break 
the force balance so that some magnetic held line tubes will rise, together with the beaded 
matter, into the magnetosphere. Note that the magnetic held of ~ 10^® G easily overcomes 
the weight of the column of 10^^ g cm“^ so that an upper layer of width ~ 100 m may be 
expelled from regions with an appropriate structure of the held. If the fraction ( of the giant 
hare energy, E = erg, is transferred to the ejected matter, a mass of lO^^C -^46 g Kiay 

be ejected. As discussed in Granot et ah (2005), if all of the inferred ejecta were released 
from the surface of the NS during the initial “hard spike” (< 0.5 s) of the giant hare, the 
outhow would be opaque to y-rays and the Dec. 27 hare would not have been observed. This 
can be avoided if there are regions on the magnetar surface from which radiation is expelled 
without matter, and other points from which matter is expelled. 

One possible observational signature of this process is the detection of ultra high-energy 
(UHE; Eu > 1 TeV) neutrinos from SGR 1806-20 coincident with the Dec. 27 hare. In 
this non-relativistic wind, internal shocks produced by signihcant variations in the outhow 
velocity within 0.5 light seconds of the star will accelerate some protons to energies high 
enough that they create pions through collisions with other protons. When these pions 
decay, they can produce TeV neutrinos. If the total energy in neutrinos is e^E, where E 
is the initial kinetic energy of the ejecta as estimated in Equations (9) and (13), then the 
observed huence of neutrinos, Eu, is: 

Eu ^ 1.2e,^^_ii744.5(i75^ X 10“^ ergs cm“^ (14) 

where ^ 44,5 = E/IO"^^-^ ergs and = e,^/0.1 (e.g. Eichler & Schramm 1978). If ~ 0.1, 
this is much higher than the 10“^ erg cm“^ typically expected from bright GRBs (Eichler 
1994). Depending on the exact values of and E, these neutrinos could possibly have been 
detected with current arrays, and the Dec. 27 event thus makes the best test case so far for 
testing the hypothesis of UHE neutrino emission from 7 -ray outbursts. It is not expected 
than any UHE neutrinos will be produced in the forward shock generated by the outflow as 
it expands into the ISM (Fan et al. 2005). 
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Table 1. 


Radio Observations at 4.8 GHz of the radio nebula produced by SGR 1806-20 


Average Epoch [UT] 

Days after Burst 

Telescope 

S' 4.8 GHz [mJy] 

2005 Jan 03.83 

6.93 

VLA 

80±1 

2005 Jan 04.61 

7.71 

VLA 

66±3 

2005 Jan 05.26 

8.36 

ATCA 

60±1 

2005 Jan 05.66 

8.76 

VLA 

57±3 

2005 Jan 05.85 

8.95 

ATCA 

53±1 

2005 Jan 06.24 

9.34 

ATCA 

46±2 

2005 Jan 06.84 

9.94 

ATCA 

39±2 

2005 Jan 06.84 

9.94 

VLA 

39±1 

2005 Jan 07.90 

11.00 

VLA 

28±2 

2005 Jan 08.19 

11.29 

ATCA 

25±2 

2005 Jan 09.07 

12.17 

ATCA 

21±1 

2005 Jan 10.07 

13.17 

ATCA 

17±1 

2005 Jan 12.06 

15.16 

ATCA 

12±1 

2005 Jan 14.08 

17.18 

ATCA 

10±1 

2005 Jan 16.08 

19.18 

ATCA 

7±1 

2005 Jan 18.01 

21.11 

ATCA 

6.5±0.5 

2005 Jan 20.01 

23.11 

ATCA 

5.5±0.5 

2005 Jan 22.08 

25.18 

ATCA 

4.5±0.5 

2005 Jan 23.08 

26.18 

ATCA 

5.5±0.5 

2005 Jan 24.62 

27.72 

VLA 

5.0±0.2 

2005 Jan 24.81 

27.91 

ATCA 

4.4±0.5 

2005 Jan 25.99 

29.09 

ATCA 

5.5±0.5 

2005 Jan 27.99 

31.09 

ATCA 

5.8±0.5 

2005 Jan 29.99 

33.09 

ATCA 

5.5±0.5 

2005 Jan 31.82 

34.92 

ATCA 

6.0±0.5 

2005 Feb 01.82 

35.92 

ATCA 

5.2±0.3 

2005 Feb 02.82 

36.92 

ATCA 

5.8±0.4 

2005 Feb 03.59 

37.69 

VLA 

4.8±0.2 

2005 Feb 03.82 

37.92 

ATCA 

4.8±0.3 

2005 Feb 05.91 

40.01 

ATCA 

4.4±0.3 

2005 Feb 07.53 

41.63 

VLA 

4.1±0.2 

2005 Feb 10.52 

44.62 

VLA 

3.9±0.2 

2005 Feb 11.92 

46.02 

ATCA 

3.6±0.4 

2005 Feb 12.62 

46.72 

VLA 

4.2±0.2 

2005 Feb 14.80 

48.90 

ATCA 

3.5±0.2 

2005 Feb 19.54 

53.64 

VLA 

3.3±0.2 

2005 Feb 20.98 

55.08 

ATCA 

2.9±0.4 

2005 Feb 21.61 

55.71 

VLA 

3.3±0.1 

2005 Feb 23.99 

58.09 

ATCA 

2.8±0.5 

2005 Feb 26.55 

60.65 

VLA 

2.7±0.1 

2005 Feb 28.85 

62.95 

ATCA 

2.5±0.4 


Note. — Flux densities before 2005 Jan 18.01 are also reported in the 
Supplementary Section of Gaensler et al. (2005b). 
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Time since flare (days) 

Fig. 1.— The 4.8 GHz light curve of the radio nebula associated with SGR 1806-20 up to 
day 62 (2005 Feb 28) after the giant flare. The circles represent data taken with the VLA, 
and stars data taken with the ATGA. The dot-dashed line in the light curve are the result 
of htting the data to the model described in §3 and whose parameters are given in the text. 
The dotted line shows the power-law component of the model £t while the dashed line shows 
the additional component due to the swept-up, shocked, ambient material. The ends of the 
dot-dashed line correspond to the hrst and last data points included in the fit. 
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Fig. 2.— Evolution of the expanding shell. The solid red line corresponds to the dimension¬ 
less radius, r, the dashed green line to the velocity v in units of the initial velocity, and the 
dotted blue line corresponds to the dimensionless synchrotron flux, /. 




